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ABSTRACT 



Some of the methods of automatic control are used to investigate 
the feasibility of a suspension system using a spring with two spring 
rates. Following a disturbance the spring rate is changed discont inuous- 
ly from one value to the other in such a manner that a damping action 
results. A single degree of freedom corresponding to the vertical 
motion of the automobile body with respect to one wheel is considered 
and this system is simulated on the analog computer. The objectives of 
this investigation are: (1) to determine the effect of the spring in 

reducing the frequency response resonance peak. (2) to determine the 
transient response characteristics of the system when subjected to a 
step disturbance. The results show that the two-rate spring does provide 
a significant amount of damping. 

The writer wishes to express his appreciation for the suggestion of 
this thesis topic and for the encouragement and a ssistance given him by 
Doctor G. J. Thaler of the Electrical Engineering Department of the 
U. S. Naval Postgraduate School. 
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DESCRIPTION 

Input signal amplitude (sine wave or step)- feet 

Spring-rate return - per cent or pounds per foot 

Viscous damping constant - pounds per foot per second 

Position error ( A measure of spi ing compression or tension) 
feet 

Rate of change of position error - feet per second 

Frequency - cycles per second 

Viscous damping force - pounds 

Spring force pounds 

Spring rate - pounds per foot 

One-fourth of the spring mass - slugs 

LaPlace Operator 

Time - seconds 

Input displacement (wheel displacement) - feet 

Output displacement (sprung mass displacement ) - feet 

Output velocity - feet per second 

Damping ratio of linear system 

Radian frequency - radians per second 

Natural frequency 

Resonant frequency 
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CHAPTER I 



INTRODUCTION 

Following the acceptance of the automobile as a reliable means of 
transportation manufacturers have continually striven to improve the 
suspension system to increase passenger safety and comfort. The obvious 
method of determining the ride characteristics of a proposed suspension 
is to install the suspension and drive and ride in the automobile over 
various types of roads to observe its performance. To determine what 
changes should be made one makes the changes and observes their influence 
on the ride characteristics. This is the observation method of ride 
analysis and this method was largely used in evolving the suspension sys- 
tem of todays automobiles. 

However with the observation method the observer can only observe 
the characteristics of one automobile at a time and therefore he must 
rely on his memory to compare one automobile suspension with the next. 

In other words the observation method is subjective and it does not furnish 
numerical data which can be evaluated and compared over long-time periods. 

Mathematical analysis using an analog computer can provide numerical 
data and the means for comparing and evaluating the characteristics of 
proposed designs. If this analysis proves a particular design to be 
feasible, then of course, the subjective observation method can be used. 
Often however, the time and expense associated with fabricating the 
suspension and conducting the observation test can be saved when the 
mathematical analysis proves the design not feasible. There is another 
advantage to the mathematical analysis. The differential equations that 
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describe the auto suspension system are of the same form as those of a 
position servomechanism. This suggests that the techniques of automatic 
control might be used in analyzing and improving the ride characteristics 
of automobile suspension systems. 

This thesis presents the results of a mathematical study of a suspen- 
sion system in which the usual constant rate spring is replaced with a 
non-linear spring; i.e., a spring with a spring rate that is switched to 
either of two values. Switching is a function of spring displacement and 
rate of change of spring displacement. Investigation is restricted to one 
degree of freedom corresponding to the vertical displacement of the auto 
body with respect to the wheels. From the system differential equations 
a block diagram is formed similar to a servomechanism and fundamental auto- 
matic control techniques are used to obtain the closed loop frequency re- 
sponse on the analog computer. Data is obtained for systems employing the 
non-linear spring with and without a viscous damping shock absorber. The 
effect of varying the viscous damping is investigated and also the effect 
of changing the ratio of the two values of spring rate. 

Chapter II presents the justifications for treating the suspension 
as a position servomechanism, and develops the control type block diagram 
that represents the automobile suspension. The defining relationships 
for the two-rate spring are given together with the reasons for consider- 
ing such a suspension spring. 

In order to compare the proposed suspension with a conventional sus- 
pension system the characteristics of a typical automobile suspension are 
explained. A summary of the literature on human comfort criteria is in- 
cluded and compared with the typical suspension system to determine how 
the suspension characteristics should be changed to improve the riding 
quality. 
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For a linear approximation to the behavior of the two-rate spring 
the describing function is derived from a Fourier analysis of the spring 
force output wave form resulting from a sinusoidal displacement function. 

To simulate the proposed suspension system on the analog computer 
is was necessary to devise a means of simulating the two-rate spring. 
Whereas analog computer literature contains methods for simulating many 
common non-linearities none of these techniques are suitable for simu- 
lating the two-rate spring. Chapter V describes the special transistor/ 
relay logic circuit designed to perform the spring rate switching opera- 
tion. 

Chapter VI and VII, present the results of frequency response and 
transient response tests respectively. Transient response data is dis- 
played in the phase portrait form. 

The final chapter summarizes the findings of the feasibility study 
and suggests further suspension system research. 
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CHAPTER II 



THE AUTOMOBILE SUSPENSION AS A SERVOMECHANISM 

The first applications of mathematical analysis to automobile sus- 
pension system performance simply consisted of solving the differential 
equations of motion and determining how stability and ride characteristics 
were related to the coefficients and form of the equations. The character- 
istics of individual system components and the constants representing 
physical parameters do not appear explicitly in the equations and there- 
fore the solutions to the differential equations of motion do not readily 
indicate how the physical system should be changed to improve the dynamic 
characteristics of the suspension. Whereas this method of mathematical 
analysis can be used to evaluate the characteristics of a particular 
suspension, trial and error techniques must be used to design system 
improvements. 

The pioneer automatic control systems engineer was faced with this 
same design handicap. The differential equations describing a control 
system do not directly reveal what changes must be made to the physical 
system to obtain a desired change in the response. For simple control 
systems trial and error methods can be used to ultimately obtain the de- 
sired response. However, as control specifications demand greater speed 
and/or greater accuracy the system must grow in size and complexity to 
include more system components, and as the system becomes more complex 
it becomes increasingly more difficult to apply trial and error techniques 
to solve the design problem. During World War II military requirements 
for high performance such as in anti-aircraft fire control systems resulted 
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in control systems that were so complex that trial and error techniques 
were no longer practical . 

To facilitate the design of these complex systems there was develop- 
ed an operational method of mathematical analysis employing La Place 
transforms* It is still a method of solving the differential equations 
of motion but an advantage results from a different approach to the solu- 
tion* Emphasis is placed on the physical arrangement of the system and 
the effect of individual components on overall system performance is more 
readily determined. The operational approach to the analysis makes possi 
ble a block diagram representation of the system that is also useful in 
relating overall system performance to the characteristics (La Place 
transform functions) of individual components. Numerous other control 
design devices are made available through the application of La Place 
transforms to the differential equations of the system. A detailed 
explanation of this operational transfer function approach and associat- 
ed techniques is available in any feedback control systems text and de- 
tails will not be given here. 

The similarity of the differential equations of motion for the auto- 
mobile suspension system to those for a position servomechanism suggests 
that the operational method and associated feedback control techniques 
might be used for analysis of suspension systems. A conventional auto 
suspension system and the equivalent block diagram developed from the 
transformed differential equations are shown in Figure 2-1. Of course, 
the complete suspension system is a five mass system (body plus four 
wheels) and if all modes of translation and rotation are considered for 
each mass the system then has 30 degrees of freedom. However, ride 
characteristics and stability are predominantly dependent upon vertical 
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displacement of the auto body and therefore this thesis work is con- 
fined to a single degree of freedom corresponding to the vertical dis- 
placement of the body with respect to one wheel. 

The differential equations for the system of Figure 2-1 are stated 
as follows. For a vertical wheel disturbance, X^(t), a spring force, 
F^Ct), result proportional to the relative displacement of the 

body with respect to the wheel. 

fm = /<[xfi)-xM H 



The viscous damping action of the shock absorber develops a damping force 

F c (c>: 



F c (t) -- -xj)] M 

The combined forces accelerate the body, M: 

F K (t) + F c (i) = MX z (i) (x-3) 



Substituting (2-1) and (2-2) in (2-3): 



K{x,ft)-m + c[xfi) -x/0}^ m x-(t) H 



To transform E quat ion (2-4) to La Place operational notation initial condi- 
tion values of X^(t) and X^(t) must be known. If the system is initially 
in vertical static equilibrium then: 



x/o] = X : (b) = X, (o) 



After transformation Equation (2-4) becomes: 
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Position error is defined as the relative displacement of the wheel with 
respect to the body. 



From the transformed equations the equivalent block diagram is 
easily developed and the similarity to a position servomechanism be- 
comes apparent. The input and output of this servo- like system are the 
vertical displacements of the wheel and body respectively. E corres- 
ponds exactly to the position error of a servomechanism. This error is 
acted on by the system gain, K/M, and the transfer function 1 + C/K S, 
until the error is reduced to zero. Note that if the shock absorber is 
removed C equals zero and the system then lacks a device for dissipating 
energy. Following a vertical disturbance imparted to the wheel, the 
energy introduced remains within the system and the body, M, will con- 
tinue to oscillate as the energy is transferred back and forth between 
the potential energy of the spring and the kinetic energy of the moving 
body mass, M. To prevent this sustained oscillation the viscous damper 
shock absorber is included in present automobile suspensions. 

In the suspension system unlike the usual servomechanism, it is not 
required or desired that the output closely follow the input. Large 
"position errors" are generally desired and slow follow up is permitted 
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so that the passenger is not subjected to displacements and accelerations 
corresponding to each and every imperfection when driving over washboard 
and pot-hole road surfaces. In other words it is desired that the output 
of the suspension system follow some general average trend of the input. 

From servo theory it is known that this low fidelity can be obtained by 
making system gain* K/M, small which means a weak spring and/or a large 
body mass. However, as the spring constant is decreased the tendency 
for the spring to "bottom 11 becomes greater, i.e., for a large disturbance 
the amplitude of oscillations exceeds the range of spring compression or 
the clearance between chassis and body. When this occurs any further 
wheel displacement is transmitted directly to the body with the obvious 
detrimental effect on passenger comfort, and perhaps structural damage 
will occur. 

Tf K/M must be small for low fidelity considerations then the tendency 
to bottom must be controlled by viscous damping, or conceivably by some 
other scheme not yet developed. The trend in automobile suspensions over 
the last two decades has in fact been toward ever decreasing K/M neces- 
sitating ever increasing shock absorber viscous damping. However, viscous 
damping introduces adverse characteristics with the useful amplitude damp- 
ing. In the transfer function K+CS it is noted that the shock absorber is 
a velocity sensitive element. Damping force, F^(s), is proportional to CSE 
where S is the La Place notation for the first derivative. Therefore, even 
small amplitude disturtances when suddenly applied cause large accelerating 
forces. 

To summarize, a suspension employing linear viscous dampers and springs 
cannot provide ride characteristics that are optimum in all respects. 
Compromises must be made and usually the least satisfactory character of 
the ride is the disturbance transmitted to the passenger through the 
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velocity sensitive shock absorber when the automobile passes over step- 
like pot-hole road imperfections. 

The insight gained by a control systems approach to suspension 
analysis suggested to Martin and Jeska^ that one means of diminishing 
passenger disturbance when riding over pot-holes etc., is to use a non- 
linear viscous damping action. The system they investigated incorporates 
a large and a small viscous damping rate. The low viscous damping rate is 
effective when the suspension is in vert ical static e 4 uilibrium and remains 
effective when a step disturbance is applied. Initial accelerating forces 
transmitted to the passenger are thus minimized. Following the step dis- 
turbance the large viscous damping rate is introduced to quickly dampen 
oscillations. 

Another suggestion was made by Professor A. G. Thompson of the Univer- 
sity of Adelaide, Adelaide, South Australia. Based on the control systems 
approach Thompson proposes a suspension system in which the viscous damper 
is omitted and oscillation damping is provided by a two-rate spring. When 
a disturbance sets the system into oscillation the spring rate is switched 
between a small value and a large value. The switching points are determin- 
ed by position error and rate of change of position error. 

Figure (2-2) illustrates the switching action. Briefly stated when- 
ever position error and rate of change of position error have like signs 
then the large spring rate is switched in and when position error and rate 
of change of position error have opposite signs the small spring rate is 
switched in. For a sinusoidal position error there will be four switching 
points each period. The large spring rate switches in as position error 



^G. J. Martin, and R. D. Jeska, IRE Convention Record, Pt I, pp 89-98, 

1953. 
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passes through zero and the small spring rate switches in as rate of 
change of position error passes through zero. Because switching occurs 
at quarter cycle intervals, Thompson calls this switching action, quarter 
cycle return. 

Thompson’s idea is based on a fundamental method of compensating a 
servomechanism to make it less oscillatory. If a network or a device is 
inserted in the forward signal path to introduce a phase advance it can be 
seen by examining the transfer function of the device that the phase advance 
accomplishes an approximate differentiation. The phase lead operation is 
effective therefore in damping out oscillations but hopefully in the auto- 
mobile suspension the approximate differentiation will not transmit high 
accelerating forces to the passenger when step inputs are applied. 

For quarter cycle spring return it is obvious that the fundamental 
component of the spring force wave form is phase advanced with respect to 
the position error signal. 

Figure (2-3) presents a more intuitive explanation of the damping 
action provided by spring rate switching. The upper panel of the figure 
shows a sketch of position error versus time for a step input disturbance. 
The lower panel schematically depicts the condition of the spring (stretch- 
ed or compressed) and which spring rate is effective at times correspond- 
ing to those labeled on the error sketch. A study of the schematic in 
conjunction with the error sketch reveals that whenever the body mass, M, 
is moving away from steady state position the large spring rate is ef- 
fective and when the body moves toward steady state the small spring is 
effective. At extremes of error all system energy is contained in the 
spring which is either in maximum compression or maximum tension. As 
maximum compression or tension is approached the large spring rate is 
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effective. Switching occurs at the maximum condition and part of the 
energy is dissipated when the small spring rate is switched in. On 
each successive swing there is less energy available to transfer to 
kinetic energy of motion. 
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CHAPTER III 



PASSENGER COMFORT CRITERIA 

The previous chapter mentions one major difference between the 
suspension system and the usual servo-mechanism. Servomechanism specifi- 
cations normally require that the output position closely follow the in- 
put position, whereas in the suspension system large position errors are 
desired to give the passenger a smooth ride over road imperfections. 

Stated another way it is desired that the verticle motion of the auto 
body only follow some average trend of the road surface profile. The 
difficulty in determining specifications for an optimum suspension system 
partly lies in the interpretation of what is to be considered the average 
trend for passenger comfort purposes. 

Measurement of the subjective effects of vertical motion is compli- 
cated by psychological effects and physical effects caused by other factors 
such as air temperature and humidity, engine noise, and other less tangi- 
ble effects that contribute to the passenger’s sense of security and well- 
being. Comfort in its broadest sense is a measure of the passenger’s free- 
dom from both mental and physical distur >ance and it is not easy to complete- 
ly separate discomfort caused by vertical motion from discomfort from other 
sources. Also, threshold of discomfort varies between men and women, with 
age, and even between two persons of the same age and sex. 

Although the experimental work of this thesis is concerned only with 
vertical motion in a system with a single degree of freedom, a complete 
study must consider passenger disturbances acting in a horizontal plane 
and the various rotational motions associated with cornering and braking. 
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Considering only vertical disturbances the ride motion can be separat 



frequency band extends from zero to about one cycle per second for a typi- 
cal automobile. In this range there is no ride smoothing,, The wheel and 
the body both follow the road surface profile. Ride quality may be im- 
proved by narrowing this frequency band but there must always be a low 
band within which no smoothing is provided so that the vehicle can follow 
the road contours of humpback bridges, large dips, and sine waves with 
crest to crest distances large with respect to the wheel base. 

The mid- frequency band extends up to about two cycles per second and 
includes the natural frequency of the body mass/spring system. In the 
closed loop transfer function of Figure (2-1), the constant term in the 
denominator determines the natural body frequency. 



For disturbances in the mid range the wheel follows the road contour and 
the body mass/spring combination amplifies the vertical disturbance and 
so the body oscillates with greater amplitude than that of the road im- 
perfections. The viscous damper shock absorber reduces this resonance 
peak but the amount of damping required to completely eliminate the peak 
will make the step response characteristic unacceptable. 



H. L. Cox, The Riding Qualities of Wheeled Vehicles, Proceedings of 
Auto Division of Inst, of Mech Engineers, 1955-56. 

3 

David Hodkin, Progress in Comfort, Part One, The Autocar, pp 747-751 
May 12, 1961. 

4 

R. Schilling and H. 0. Fuchs, Modern Passenger-Car Ride Character- 
istics, J. Appl . Mech., 8, pp A-59-65, June, 1941. 
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The body is essentially unaffected by road disturbances in the upper 
frequency band above two cycles per second. In this band the filtering 
action of the body mass/spring combination attenuates the amplitude of 
wheel oscillations. The wheel continues to follow the road contour 
when the disturbing frequency falls in the lower end of this band. How- 
ever, for disturbances with a frequency near the natural frequency of 
the unsprung mass between the tire and the suspension the amplitude of 
wheel oscillation exceeds the amplitude of road surface contour. This 
natural wheel resonance frequency is a function of the unsprung mass and the 
spiihg rate of the tire and equals about eight cycles per second for a 
typical automobile. The viscous damper shock absorber reduces this wheel 
resonance but even without damping, wheel resonance is normally sufficient- 
ly beyond the pass band of the suspension system so that passengers do not 
feel wheel resonance vibrations. (Also wheel resonance normally occurs 
below the audio threshold.) For still higher frequencies ride smoothing 
action is accomplished entirely by the tires. The tire follows the road 
imperfection while the wheel and body ride smoothly. 

In the above discussion we have considered the automobile suspension 
to be excited by a single frequency. That is ; constant speed over uniform- 
ly spaced ripples has been assumed. The frequency then is proportional to 
the speed and inversely proportional to the crest to crest separation or 
wave length. Occasionally a uniformly spaced ripple or washboard road 
surface is encountered. However road surface imperfections are usually 
randomly spaced with varying amplitudes, and consequently the suspension 
is simultaneously subjected to disturbances falling in all three frequency 
bands. 

With the fundamental suspension characteristics established qualita- 
tively in terms of a frequency spectrum, physiological response to 
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sinusoidal motion will now be considered. A search of the literature 
revealed numerous articles on the subject of human comfort criteria for 
harmonic motion. Although at first glance one might expect a comfort 
criteria for sudden displacements such as step functions* the complete 
"passenger suspersion" system includes the seat cushions in addition to 
the auto suspension. Therefore to subject the passenger to a vertical 
step input, the automobile should be dropped from a tall building or be 
subjected to some other unlikely disturbance. Two of the more complete 
quantitative treatments of vertical vibration human comfort criteria will 
be summarized. 

Lippert, a research analyst for Douglas Aircraft Company, has recon- 
ciled graphically the findings of a number of studies on the effect of 
vertical vibration on human comfort.^ Figure (3-2) presents his daca 
for thresholds of perception and discomfort. It is likely a common 
experience that the comfort reaction is predominately related to accelera- 
tion and jerk, i.e., to the second and third derivative of displacement. 
For simple harmonic motion the displacement is of the form: 



X-As/m/jtt (?-:) 

The first, second, and third derivatives are respectively 

velocity = X = ZAfft C.'\T 1/7‘7 t ( 3 - r) 

acceleration = X — "'''/ '\! 7 2 / Si/'( ,17 1l ' (~A J 

jerk = X — ' SAV 1 '~ f COSZTtft (3- - S) 



Because acceleration and jerk are proportional to the second and third 
power of frequency respectively by intuition we would expect threshold 
amplitude to decrease with increasing frequency as found on Figure (3-2). 

^Lippert, Human Response to Vertical Vibration, SAE Journal, 55, 
pp 32-34, May 1947. 
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Janeway, Head, Dynamics Research Department, Chrysler Corporation 

has made a similar study and presents his data in a slightly different 

6 

manner. The vibration spectrum is divided into three frequency bands 
in each of which a different measureable characteristic of vibration 
determines the threshold level* 

The low range extends from 1 to 6 cycles per second. In this range 

jerk is the limiting characteristic of vibration and the comfort limit 

3 

is established at 40 Ft/Sec„ An equivalent limit expression is: 

Af — 2 • (A~ <amp//fuc/e m iak/icsj (3-6/ 

( f - freouehcy //, cps) 

In the mid-range extending from 6 to 20 cps acceleration is the 

limiting characteristic and the threshold of discomfort is established 

2 

at 3.37o g - 1.07 Ft/Sec. or: 

Af' ^ 3 (3-j 

Velocity determines the threshold in the high frequency range, 20 
to 60 cps. Maximum velocity amplitude is .105 inches per second. 

Af — 2o (p-e) 

A comparison of Lippert*s data with Janeway p s shows that they are 
reasonably consistent. Lippert f s piecewise linear representation places 
the threshold of discomfort at slightly lower amplitude values than does 
Janeway’s closed form relationship, in Equations (3-6, 7,8). 



R. N. Janeway, Passenger Vibration Limits, SAE Journal, 56, pp 48- 
49, August 1948. 
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Comfort criteria can \ ow be correlated with suspension system frequency 
characteristics of Figure (3-1) . Body resonance falls within the jerk 
sensitive range of comfort. To improve the frequency response in this 
range, it is desired to reduce the resonance peak and lower the resonant 
frequency. 

Wheel resonance falls within the acceleration sensitive range. Here 
also it appears desirable to lower the resonant frequency to take advan- 
tage of a higher discomfort threshold. However wheel resonance must be 
kept beyond the suspension band-pass cutoff. 

The comfort criteria presented in this chapter are used in the 
experimental work as a means of comparing the two-rate spring suspension 
with the normal suspension. 
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CHAPTER IV 



DESCRIBING FUNCTION FOR THE TWO- RATE SPRING 

The suspension system investigated has been made intentionally 
non-linear by including a spring with a spring rate discont inuously 
switched between a larger and a smaller value. This makes the dif- 
ferential equations of the system non-linear. Likewise the transfer 
function of the spring is non-linear. Usually non-linearities make the 
differential equations unsolvable by classical mathematical methods and 
automatic control methods used on linear systems are often awkward or 

not applicable. This chapter treats the non-linear spring by means of 

7 

describing function theory, which is a control systems method developed 
to extend linear transfer function techniques and frequency response 
techniques to systems containing a non-linear component. 

Brief 1} stated a describing function is an approximate linear trans- 
fer function used to represent a non-linear component. If the approxi- 
mation is a good one then all of the linear system frequency response 
techniques can be used to analyze the non-linear system. The describing 
function relates the fundamental component of output waveform of the non- 
linearity to the input. The following assumptions are made. (1) the in- 
put waveform is a pure sine wave; (2) the output waveform can be approxi- 
mated satisfactorily by the fundamental component of the Fourier series 
representation; (3) no zero frequency component or sub-harmonics are 
present . 

{ 2 . J. Thaler and M. P. Pastel. Analysis and Design of Non-linear 
Control Systems. McGraw-Hill Book Company, Inc., New York, 1962. 
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For the two-rate spring non-linearity the spring force output waveform 
is conveniently expressed by a Fourier series. Therefore, the describing 
function is obtained simply by taking the quotient of the fundamental 
term of the Fourier series and the input sine wave. To write the Fourier 
series for the output waveform it is helpful to express the defining re- 
lationships of Figure (2-2) in another form: 



K zz. tahgeh sph/flg /<?£> 
f(~3 — Stmi/zt spt-wg No ? • O -^b b 

Then the defining equations are: 

F (s) “ /< t (s) (uh th 

F(s) — (K~B) ECs) t*)hen 




Here B is a measure of the amount of spring rate return, i.e., the 
amount by which the larger rate is reduced by switching . Using Equations 
(4-1) and (4-2) to obtain the output Fourier series the describing func- 
tion, G^, is found to bet* 




*The describing function is derived in Appendix I 
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In general, describing functions may be used to predict the closed 
loop frequency response and to investigate system stability. In this 
thesis work the describing function is used only as a preview to the 
behavior of the two-rate spring. Equation (4-3) confirms the suspicion 
that the two-rate spring transfer function is insensitive to input amplitude 
and frequency. For the frequency response experimental procedure this 
means that the input signal to the spring transfer function, G^, need not 
be held constant as the input frequency is varied to produce the frequency 
response data. Equation (4-4) indicates that the two-rate spring intro- 
duces a constant phase advance independent of frequency. Recall that it 
is Thompson's desire to introduce such a phase shift to create an approxi- 
mate derivative effect* Equation (4-4) also shows that the linear approxi- 
mate transfer function produces the maximum phase advance for 100% spring 
return, i.e., for the limiting case B equals K for a maximum phase ad- 
vance of arctan 2/n = 32.6 degrees. 
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CHAPTER V 



ANALOG SIMULATION OF THE TWO- RATE SPRING 

To simulate the proposed two-rate spring suspension on the analog 
computer it was first necessary to find a technique for simulating the 
switch action of the two-rate spring. In the language of the computer 
the gain of an operational amplifier must be increased and decreased dis- 
continuously as prescribed by the defining equations of the spring. Where- 
as analog computer texts and the literature describe techniques for simulat- 
ing the common non-linearities such as saturation, relay dead-zone, relay 
hysteresis etc., none of these techniques are suitable for the two-rate 
spring. This chapter describes the circuit designed to perform the switch- 
ing operation. 

A written statement of the switching required suggests the use of 
logic circuits. 

When position error and rate of change of position error 
are both positive or when position error and rate of change 
of position error are both negative then the larger spring- 
rate is effective. 

Two " and " logic circuits and one "or" logic circuit are required to 
determine switching conditions. Amplifier gain can be easily changed by 
switching in a resistor. Therefore, a parallel combination of two relay 
contacts can be used for the or logic. If normally open relay contacts 
are used any of the four arrangements in Figure (5-1) can be used. The 
particular " and " logic circuits chosen require that either (4-1-b) or 
(4-1-c) be used to provide the larger amplifier gain when the "or" condi- 
tions are fulfilled. 

Diode " and " logic circuits were first tried for operating the M or" 
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relays. However the diode ’’and” circuit output current level is the 
same as the input signal current level which in this case is only the 
current available from an operational amplifier. Pull-in current for 
the available relays exceeds the output of an operational amplifier and 
therefore the diode logic was abandoned in favor of transistor logic. 

A schematic of the transistor "and" logic circuits is shown in 
Figure (5-2). Each relay coil is operated by a two transistor "and" 
circuit. PNP transistors are used in a common emitter configuration. Both 
collectors are tied together and connected to one end of the relay coil. 

Each "and" circuit operates as follows. A positive input to the base of a 
transistor turns that transistor off and a high resistance appears in the 
emitter - collector path. A negative base input turns that transistor on 
and emitter-collector resistance drops to several ohms effectively clamping 
both collectors to giound. Therefore both transistors must be turned off 
with positive inputs to energize the relay. As a result when E and E are 
both positive, relay number one closes, and when -E and -E are bot\ posi- 
tive then relay number two closes. In this way spring-rate switching is 
accomplished as called for by the defining equations. 

Some refinement is necessary to allow correct relay switching indepen- 
dent of position error amplitude. (Position error for a step input is a 
decaying sine wave oscillation). It is desired that each transistor turn 
fully on, i.e., saturate, whenever the base signal passes through zero 
going negative and conversely each transistor should turn off immediately 
when the base signal passes through zero going positive. Therefore it is 
necessary to use an amplifier with high gain and diode limiting to pro- 
duce an approximately square position error wave. Likewise rate of change 
of position error is produced with a diode -limited high gain differentiating 



28 



AN D —Ecc o * AN D 



WH" 



V. 'V 




QC ^ 




£ X 

8 
X 



£ 



d 2 
<$: ^ 





amplifier. These auxiliary amplifiers are shown in Figure (5-3), The 
resistor in series with the input capacitor of amplifier A-3 depresses the 
noise pick up inherent with the differentiating circuit. Unfortunately the 
series resistor also introduces an undesired frequency and amplitude sensi- 
tive phase shift in the square wave output. For correct switching through- 
out a frequency response test it is necessary to adjust the capacitance value 
to compensate for phase shift each time the frequency is changed. 

As discussed in a later chapter, midway through the step response 
tests it was found preferable to redefine the switching requirement in terms 
of rate of change of output position instead of rate of change of position 
error. Thus the dif ferentiat itn amplifier is eliminated and the phase shift 
problem is circumvented. 
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CHAPTER VI 



FREQUENCY RESPONSE TESTS 

To have a basis for qualitatively comparing the two- rate spring 
suspension with a typical conventional suspension (Figure 3-1), closed 
loop frequency response tests were performed* The block diagram of the 
two-rate spring suspension system (Figure 6-1) is identical with that for 
the conventional suspension (Figure 2-1) except that the spring must now 
be represented with the describing function, G^, or by the defining equa- 
tions (4-1) and (4-2) instead of the constant numeric K. The larger spring- 
rate is chosen as 720 pounds per foot for all tests and M, one fourth of 
the sprung mass is 18 slugs. (These values are typical of a light weight 
passenger automobile although for this qualitative investigation numerical 
values are not critical. They do not change the character of the response 
but only serve as a scaling factor and here the interest is in a comparison 
of the responses with and without spring rate switching.) Viscous damping, 

C, is one of the variable parameters and for each test C is chosen to give 
a desired damping ratio. Throughout the remainder of the thesis "damping 
ratio 11 implies the damping ratio of the linear system without spring-rate 
switching. 

Dcu-np/hg P\abo — ^ — ^ZfMK (g~i) 

B, spring-rate return^ is the other variable parameter. Choosing B determines 
the value of the smaller spring rate since the larger spring rate is con- 
stant at 720 pounds per foot. 

The Donner Model 3400 analog computer was used and Figure (6-2) shows 
the problem board circuit airangement. 
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To verify that the two-rate spring characteristic is not a function 
of input signal amplitude a preliminary test was made holding input fre- 
quency constant but varying input sine wave amplitude. This data, tabu- 
lated below, does indicate the spring characteristic to be independent of 
input signal amplitude. Therefore it is not necessary to hold position 
error constant throughout each frequency response test. 
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For the fast frequency response test no viscous damping was used 
(C = 0) and spring-rate return was set at B = 360 pounds per foot. This 
corresponds to 50 per cent spring return and zero damping ratio. The data 
are plotted on Figure (6-4.) 

Without spring-rate return the amplitude of the resonance peak would t 
infinite. The finite resonance peak of Figure (6-4) demonstrates the 
damping effect at the spring rate switching. The natural frequency for 
the system is 



For a lowly damped system the resonant frequency is very nearly equal to 
the natural frequency and therefore if spring rate switching were not appli 
ed to the system of Figure <6-4) then the resonance peak would occur at 
about 6.3 radians per second. Figure (6-4) indicates that with spring-rate 
switching a significant shift of the resonant frequency to the left is 
introduced. Recall that the comfort criteria presented in Chapter III 
dictates a reduction in the peak and a shift to the left to make resonance 
less perceptible to the passenger. Spring rate switching exhibits the 
ability to do both. However the resonance peak appears excessive for 50% 
spring return without viscous damping. Applying equation (3-6), for an 
arbitrary disturbance amplitude of two inches: 







<£. ■ 






The jerk discomfort threshold limit is: 
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To reduce the resonance peak, greater spring return can be used or 



some viscous damping can be added. For the next series of frequency 
response tests, spring return ; B ; was held at 360 pounds per foot while 
tests were made with various amounts of viscous damping. Figures (6-5), 
(6-6), (6-7) and (6-8) display the frequency response data for damping 
ratios of 0,2, 0,3, 0,5, 0,7 respectively. Also plotted on each figure is 
the theoretical frequency response curve for the system without spring-rate 
switching. For each damping ratio value the spring-return resonance peak 
is smaller than the corresponding linear system resonance peak. This in- 
dicates that the spring- return contributes a significant amount of the 
total damping. Spring-return resonance peak is approximately two- thirds 
as great as the linear resonance peak for each case. For 50 per cent spring 
return, to reduce the maximum jerk below the limit specified by equation 
(3-6) a damping ratio larger than 0.2 is required, based on an arbitrary 
two inch disturbance amplitude. 

Before proceeding with further frequency response testing, transient 
tests were performed. This data is presented in Chapter VII. During the 
transient testing it was discovered that spring- rate switching should be 
altered to be a function of position error and rate of change of output 
position, i.e., the modified defining equations are: 









Several frequency response tests were performed with the modified switch- 
ing. Figure (6-9) compares the original switching with the modified 
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switching for 50 per cent spring return and 0.3 damping ratio. With 
modified switching a smaller re: onance peak and a greater resonant fre- 
quency shift is evident * 

For the final frequency response tests modified switching was incorpor- 
ated with greater spring-return. Figure (6-10) presents the results for 
75% spring-return and 0.3 damping ratio, and Figure (6-11) gives the re- 
sponse for 90% spring-return and 0.1 damping ratio. As expected greater 
spring return provides greater damping and shifts the resonant peak fur- 
ther to the left. Considering only frequency response data the character- 
istic of Figure (6-11) is better than for all other parameter combinations 
investigated. Adequate damping is provided to meet the jerk comfort crite- 
ria with only a small amount of viscous damping. 
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CHAPTER VII 



TRANSIENT RESPONSE TESTS 

Servomechanisms are larely if ever excited by sinusoidal signals. 
Therefore frequency response characteristics are useful to the control 
engineer only as a means of predicting the time response of a servomechan- 
ism subjected to step* ramp s or other input signals. When possible the 
time or transient response characteristics are determined directly by ap- 
plying an input such as a step and observing the output quantity as a 
function of time. So it is with the automobile suspension system. Wash- 
board surfaces do approximate a sinusoidal disturbance but the automobile 
must ride over aperiodic road imperfections as well. Note that the fre- 
quency response curve of Figure (6-8) where a large damping ratio is em- 
ployed does not disclose the adverse derivative characteristic associated 
with suddenly applied disturbances and discussed in Chapter II. Recall 
that it is this characteristic of transmitting high accelerating forces 
to the passenger when a wheel hits a pot-hole that precludes the design 
of a linear suspension system that is optimum in all ride characteristics. 

To investigate the performance of the two-rate spring when the auto- 
mobile passes over pot-holes or rides over curbs* transient response tests 
were performed. The analog computer circuit shown in Figure (6-2) is also 
satisfactory for step response testing. However the logic switching circuit 
does not perform as well for transient testing. Whereas the differentiat- 
ing capacitor (Figure 5-3) can be adjusted for correct switching before a 
frequency run and switching will remain correct* for a step response run 
position error is a damped sine wave and the decreasing amplitude advances 
the switching time slightly as the run progresses. To minimize the effects 
of switch time advance the capacitor can be adjusted to cause slightly late 
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switching at the start of the response to compensate for advance during 
the run, 

A common method of displaying the transient response characteristics 
is on plots of output position vs. time, position error vs. time, out- 
put velocity vs. time, etc. However, the effect of spring-rate switching 
on the response of the suspension to a step input can better be analyzed 
on a plot of position error vs. rate of change of position error. Such a 
plot is known as a phase trajectory and a family of several phase trajector- 
ies, each for a different step magnitude, forms a phase portrait. Inter- 
pretation of the experimental data is enhanced by a preliminary discussion 
of the expected shape of the phase trajectory. 

From the initial value theorem the values of position error and rate of 
change of position error immediately following application of a step are:* 




On the sketch of Figure (7-1) the above initial values are the coordinates 
of point (0) „ The trajectory starts at the origin before the step input 
and arrives at point (0) immediately after ; the path in between is undefin- 
ed. For the segment between (0) and (1) the small spring-rate is effective 
and the body swings back toward equilibrium. Note the gradual curvature 
indicating low acceleration. Passing through (1) the position error re- 
verses sign and the large spring-rate is switched in as the body overshoots 
the equilibrium position. The greater curvature in the third quadrant in- 
dicates greater deceleration resulting from the large spring-rate. E 

*The initial value calculation appears in Appendix II 
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changes sign as the trajectory enters quadrant II and the small spring 
rate becomes effective while the body moves toward equilibrium position. 

For the first series of step response tests, switching is determined 
by Thompson's original defining relations (Figure 2-2). The phase por- 
traits of Figure (7-2) and (7-3) are for damping ratios of 0.3 and 0.5 
respectively. 50 per cent spring return is used for both. An analysis 
of Figures (7-2) and (7-3) follows; 

1. For equal step inputs the larger damping ratio of (7-3) 
produces a greater initial velocity as predicted by Equa- 
tion (7-2). 

2. Curvature in quadrants IV and III is greater for (7-3) demon- 
strating larger decelerating damper force. Although initial 
velocity is greater in (7-3) the velocities of (7-2) and 
(7-3) have equalized at the end of the fourth quadrant seg- 
ment due to the larger damping ratio in (7-3). 

3. Overshoot is less for (7-3). 

None of the trajectories of Figures (7-2) and (7-3) display a dis- 
continuous slope at switching points. The next step response test is 
conducted with more spring return to make switching discontinuities apparent. 
Figure (7-4) presents the phase portrait for 75 per cent spring return and 
0.3 damping ratio. The effect of the increased spring return can be seen 
by comparing with Figure (7-2), the phase portrait for the same damping 
ratio but for 50 per cent return. 

1. For equal step inputs the initial velocities are equal. This 
is consistent with Equation (7-2) which is a function only <f 
viscous damping and the sprung mass. 
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2. In the fourth quadrant the smaller spring-rate is ef- 
fective for each portrait; 360 lbs/Ft and 180 Ibs/Ft for 
(7-2) and (7-4) respectively. For equal position errors 

(spring compressions) the smaller spring rate develops a 
smaller decelerating force and hence the fourth quadrant 
curvature is less for (7-4). 

3. In the third quadrant both systems have equal (large) spring- 
rates effective and therefore both families have the same 
curvature in this quadrant. 

4. The second quadrant trajectories of Figure (7-4) exhibit mark- 
ed discontinuities of slope. It is noted with concern that 
this discontinuity does not fall on the negative E coordinate 
axis as called for by the defining equations for spring-rate 
switching . 

An investigation of the output waveforms of all computer amplifiers 
in the circuit of Figure (6-2) explains the late switching encountered 
above. These waveforms are shown in Figure (7-5). Harmonic distortion 
introduced in the spring force wave form by switching is not well filtered 
by the succeeding circuit stages. Therefore distorted output position sub- 
tracts from the input sine wave to produce a distorted position error. When 
the slope of this distorted position error becomes zero, the logic circuit 
performs the spring-rate switching. However, because of the distortion and 
the phase shift introduced by the spring-rate return the output velocity is 
past zero when switching is accomplished. 

The above explanation of the switching difficulty also suggests a 
solution: Replace rate of change of position error with output velocity 

for logic circuit switch determination. This solution conveniently elimi- 
a£es the differentiating amplifier and with it the necessity to continual- 
ly adjust for correct switching. In the remainder of the thesis switching 
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-based on Thompsons proposal will be referred to as original switching. 
Switching as a function of output velocity and position error will be 
referred to as modified switching,, 

A final series of phase portraits were obtained using modified switch- 
ing. Figures (7-6) and (7-7) present the results ‘for damping ratios of 

\ 

V s * 

0.3 and 0.5 both with 507o spring- return. The. increased damping effect 

obtained with modified switching is apparent in (7-6) and (7-7). A 

* 

physical explanation can be given for the reduced damping that results 
from original switching. As the trajectory approaches the negative E 
coordinate axis the large spring-rate is effective. Maximum compression 
is attained at cross-over^ and because the velocity instantaneously is 
zero all of the system energy is contained in the compressed spring. If 
the spring rate is switched at this point some of this energy is dissipat- 
ed and spring force decreases instantaneously to the small spring-rate 
value. If switching is late the large spring force associated with the . 
fully compressed large spring-rate is available to accelerate the sprung 
mass. Higher velocity results and when switching finally occurs much of the 
potential energy of the spring has already been transferred to the sprung 
mass and therefore a smaller amount is dissipated by the switching. Fig- 
ures (7-8) and (7-9) show phase trajectories for 0.1 and 0,3 damping ratios, 
both with 75 per cent spring-return. Figure (7-8) exhibits a desirable low 
initial velocity but overshoot is excessive and more return damping is re- 
quired. Figures (7-10) and (7-11) present the final tests. Damping ratios 
are 0.05 and 0.1 respectively and spring-rate return is 90 per cent. Both 
desirable characteristics are exhibited. Initial velocity is low because of 
small viscous damping but spring-rate return danping is adequate to limit 
the overshoot. 
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CHAPTER VIII 



SUMMARY OF FEASIBILITY STUDY 

The experimental work accomplished thus far does not constitute a 
complete feasibility study of the two-rate spring* Further frequency re- 
sponse tests and step response tests should be made and this work should 
be extended to include the effects of other components that contribute 
to "passenger suspension" such as tires*, and seat cushions* The results 
of the frequency response tests and step response tests reported in the 
previous chapters are encouraging*, 

Frequency response and transient response characteristics both demon 
strate the significant damping action introduced by spring rate switching 
Although spring-rate return appears to provide insufficient damping to 
permit the complete removal of viscous dampers at least viscous damping 
can be reduced to make step like road imperfections less perceptible to 
the passenger* The two rate spring also improves ride quality by shift- 
ing the automobile body resonant frequency toward the less perceptible 
end of the comfort spectrum* 

The following work is suggested: 

1* Conduct further frequency response tests for small damping 
ratios and large spring rate return using the modified 
switching criterion* 

2 a Extend this work to the more complete and more complex 
system including the unsprung wheel mass and spring-rate 
of the tire and perhaps the spring rate and damping of the 
seat cushion* This would involve at least two degrees of 



freedom and would permit investigation of wheel resonance 
(See Figure 3-1) „ 

3o Determine the response of the two-rate spring suspension to 
ramp inputs and combined inputs such as a positive step 
followed by a negative step* a negative ramp followed by a 
positive ramp* etc, Statistical inputs could also be consider- 
ed, 

4, Include a non-linear damper with the two-rate spring such as 
proposed by Martin and Jeska,* The attempt here would be to 
completely eliminate the initial value velocity associated 
with a step input, (See Equation 7-2) , 

5, Devise a mechanical method of accomplishing spring-rate 
switching on a real spring. 



hbid. 
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